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MODEL BASED FAULT DETECTION AND 
DIAGNOSIS METHODOLOGY FOR HVAC SUBSYSTEMS 

Abstract 

A fault detection system for an HVAC system including sensors for 
measuring the performance of a condenser, a compressor, an evaporator, and a chiller. 
The fault detection system is provided with a thermodynamic preprocessor for 
calculating characteristic quantities (CQ's) from a plurality of measured inputs. Also 
provided is a base-case lookup table for storing plural first sets of CQ values 
generated by the thermodynamic preprocessor in an initial period in which fault-free 
operation of the HVAC system is assumed, a first set of CQ values being generated 
for each of plural different measured input values. An interpolator is provided for 
interpolating a set of base-case CQ values from the first sets of CQ values stored in 
the base-case lookup table for a given set of measured inputs. A fault-detector detects 
a fault when a difference between actual CQ values and base-case CQ values exceeds 
a predetermined threshold value for at least one of the CQ values. A fault classifier 
classifies a detected fault based on which ones of the actual CQ values exceed the 
interpolated CQ values. 
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MODEL BASED FAULT DETECTION AND 
DIAGNOSIS METHODOLOGY FOR HVAC SUBSYSTEMS 



FIELD OF THE INVENTION 

The present invention relates to a method for detecting and diagnosing faults 
in HVAC subsystems. More specifically, the method of the present invention detects and 
diagnoses faults using a combination of a physical (thermodynamic) model and a neural 
network. 

BACKGROUND OF THE INVENTION 

Vapor-compression systems account for a very significant portion of energy 
consumption in the industrial and commercial sectors. In large office buildings for example, 
it is estimated that 10% to 25% of the total electric consumption can be attributed to cooling 
systems alone. Furthermore, these percentages can be much greater if a chiller subsystem 
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1 is operating at low performance levels due to the presence of faults. The extensive research 

2 in Fault Detection and Diagnosis (FDD) in HVAC systems has been motivated by concerns 

3 that range from: the need to reduce power consumption and consequently energy costs; to 

4 improve comfort levels in buildings; to reduce wear on various HVAC equipment; to reduce 

5 the magnitude of greenhouse emissions; and the need to support building operators in 

6 decision-making for building optimization. 

7 Conventional Energy Management and Control Systems (EMCS) are limited 

8 in that they do not provide the operator with the tools to diagnose faults in real-time. If 

9 effective FDD tools are designed and implemented in EMC systems, then the detection and 

10 diagnosis of various types of faults can be done automatically. 

1 1 A fault, in the context of HVAC applications, is defined as an unsatisfactory 

12 or unacceptable condition in the operation of a system or subsystem. A condition is 

1 3 unacceptable if it is a failure or if it causes one directly or indirectly through a series of other 

14 faults. There are various types of faults, some of which are more difficult to detect and 

15 diagnose than others. The three basic fault categories, ranked in order of severity, are 

16 degradations, malfunctions and hard failures. 

17 Degradation faults occur at a gradual rate but progressively worsen over a 

1 8 period of time. Overall system performance is not drastically changed until the degradation 

1 9 has matured beyond a critical level. Examples of system performance degradation faults are 

20 the fouling of condenser or evaporator tubes, or the clogging of filters. Notably, dirt and 

21 grime builds up on the heat transfer surfaces of the heat exchangers over time, which in turn 
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1 modifies the heat transfer coefficient of these devices. These faults gradually result in 

2 pressure drops, reduced flow rates, or higher temperature differences between the two fluids. 

3 Malfunction faults have more immediately noticeable end effects than 

4 degradation faults, and tend to justify more immediate service. An example of this category 

5 of faults is valve leakage, sensor errors, controller breakdowns and damper stoppages are 

6 other examples. These faults usually result in failure to maintain the desired set points in 

7 temperature, pressure or flow rates. 

8 Seized compressors and broken fan belts are examples of "hard" or "complete" 

9 failure faults. In such circumstances the impact of the fault is severe, and justifies the most 

10 immediate service. Building operators must often shut down entire subsystems in order to 

1 1 carry out the necessary repair. 

12 Fault Detection and Diagnosis (FDD) is the initial step towards taking 

1 3 corrective or preventative measures in an H VAC system. Detection and diagnosis sometimes 

14 overlap and other times are treated separately. Fault detection involves the determination 

15 that a fault truly exists based on an observable quantity exceeding some predetermined 

16 threshold or criterion. The choice of this criterion or threshold is a tradeoff between the 

17 sensitivity of the FDD scheme and the likelihood of sounding a false alarm. Fault diagnosis 

18 is the subsequent step, that isolates the cause of the fault. Sometimes this diagnostic step 

19 does not immediately locate the root cause of the fault initially detected but may involve a 

20 series of steps that eventually converges on the cause at some later point in the diagnosis. 
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1 Once the fault has been diagnosed, fault evaluation is performed to assess 

2 whether the impact of the fault sufficient to require immediate service. In general, service 

3 should be performed whenever: (1) Comfort cannot be maintained; (2) Equipment or 

v 

4 personal safety is compromised; (3) Environmental damage is done (e.g. refrigerant leakage); 

5 and (4) Service expense is justified by reduced energy costs. 

6 There are several approaches to fault detection. Most, if not all, involve 

7 observation of the differences between an "actual" quantity of the system (or component) 

8 during its operation and some "predicted" value of this same quantity. The quantity under 

9 scrutiny is sometimes the "raw" measured variables (e.g., temperatures, pressures or humidity 

10 ratios that have not been pre-processed) or it may be the result of some preprocessing that 

1 1 transforms the measured variables collectively (in specific subsets) or individually. Some 

12 examples of the quantities used for comparison in fault detection arc: Measured equipment 

13 performance vs. Model-based prediction of equipment performance; Measured equipment 

14 performance vs. Common sense expectations of acceptable performance; Measured 

15 thermodynamic states (e.g., temperatures, pressures, humidity ratios, etc.) vs. Model-based 

16 predicted thermodynamic states; Calculated (estimated) physical parameters for "current" 

17 performance vs. Calculated (estimated) physical parameters for "baseline" performance. 

18 Similarly, there are different approaches to fault diagnosis. The result of the 

19 diagnosis is not binary as in fault detection (i.e., fault or no-fault) but involves a selection 

20 from a range of different possibilities. One such fault diagnosis approach is a rule-based 

21 diagnosis, which involves a set of rules such as directional change in temperature deviations 
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1 that may be unique to a given fault. Another approach involves comparison of physical 

2 parameters for "current" operation and "normal" operation. For example, the conductance- 

3 area product, UA of a condenser can be calculated from entering and leaving temperatures, 

4 and may be used to diagnose condenser fouling. Yet another approach uses a pattern 

5 recognition technique applied to "current" residuals and a matrix of "expected" residual 

6 changes associated with each fault. 

7 Different kinds of models have been used in FDD methods. However, 

8 conventional models are unsatisfactory as they have limited general applicability. Notably, 

9 conventional models require extensive, time consuming tuning in order to tailor the model 

10 to a specific implementation. As a result, such models are not readily portable to different 

1 1 HVAC implementations. 

12 Accordingly, one object of the present invention is to provide an improved 

1 3 model-based FDD methodology that uses a combination of a physical model and a black-box 

14 model which has general applicability to HVAC systems. 

15 A further object of the present invention in this research is to develop an 

1 6 improved model-based FDD methodology for determining when and where problems occur 

17 in a centrifugal chiller system, and that is suitable for online implementation using day-to- 

1 8 day operating field data. 

1 9 Another object is to provide an improved model which compares characteristic 

20 quantities for everyday operating conditions to a predetermined baseline of a chiller' s normal 
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1 (fault- free) operation, and uses statistical hypotheses to determine if a significant deviation 

2 from this normal operation has occurred. 

3 Yet another object is to provide a FDD method which accounts for sensor 

4 errors to provide a more robust FDD methodology. 

5 SUMMARY OF THE INVENTION 

6 The above-listed objects are met or exceeded by the present fault detection and 

7 diagnosis system which includes sensors for measuring the performance of a condenser, a 

8 compressor, an evaporator, which are all components of the chiller subsystem. The fault 

9 detection system further includes a thermodynamic preprocessor, a base-case lookup table, 

10 an interpolator, a fault-detector and a fault classifier. 

11 The thermodynamic preprocessor, which is preferably a software routine 

1 2 executed on a computer or the like, calculates characteristic quantities (CQ's) from a plurality 

13 of measured inputs. The base-case lookup table stores plural first sets of CQ values 

14 generated by the thermodynamic preprocessor in an initial period in which fault- free 

1 5 operation of the H VAC system is assumed, a first set of CQ values being generated for each 

1 6 of plural different measured input values . 

1 7 The interpolator generates a set of interpolates base-case CQ values for a given 

1 8 set of measured inputs from among the first sets of CQ values stored in the base-case lookup 

19 table that is stored in memory of the computer. The fault-detector, which is preferably a 

20 software routine executed on a computer or the like, detects a fault when a difference 

6 
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1 between actual CQ values, calculated by the thermodynamic preprocessor for a given set of 

2 measured inputs, and the interpolated set of base-case CQ values for the given set of 

3 measured inputs exceeds a predetermined threshold value for at least one of the CQ values. 

4 The fault classifier, which is preferably a software routine executed on a computer or the 

5 like, classifies a detected fault based on which ones of the actual CQ values exceed the 

6 interpolated CQ values. 

7 BRIEF DESCRIPTION OF THE DRAWINGS 

8 FIGURE. 1 is block diagram showing the overall structure of the model-based 

9 fault detection and diagnosis system of the present invention; 

10 FIG. 2 is a block diagram showing the inputs and some of the outputs of the 

1 1 thermodynamic preprocessor of the present invention; 

12 FIG. 3 is an EES (Engineering Equation Solver) procedure used to determine 

1 3 the logarithmic mean temperature difference between the evaporator saturation temperature 

14 (T_evap), the chilled water supply temperature (T_chws), and the chilled water return 

15 temperature (T_chwr); 

16 FIG. 4 is an EES procedure used to determine the logarithmic mean 

17 temperature difference between the condenser saturation temperature (T_cond), the 

18 condenser water supply temperature (T_cws) and the condenser water return temperature 

19 (T_cwr) 
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1 FIG. 5 depicts equation sets for the condenser side input data used by the EES 

2 procedure of FIG. 4; 

3 FIG. 6 are EES program statements for inputting data used in the equations of 

4 FIGs. 3 and 8; 

5 FIG. 7 are EES program statements for inputting data used in the equations of 

6 FIG. 9 and FIG. 10; 

7 FIG. 8 are equations used to determine evaporator related characteristic 

8 quantities; 

9 FIG. 9 are equations used to determine compressor related characteristic 

10 quantities; 

11 FIG. 10 are equations used to determine condenser related characteristic 

12 quantities; 

13 FIG. 1 1 is an expansion valve equation used by the EES code in FIG. 8; 

14 FIG. 12 is a list of thermodynamic properties and conversion constants used 

15 by the EES code in FIGs. 8-10; 

16 FIG. 13 is a schematic of a chiller showing the location of each of the sensors 

17 used to obtain the measured inputs; 

18 FIG. 14 illustrates the structure of a single neuron unit; 

1 9 FIG. 1 5 shows the general architecture of a interconnected neural network with 

20 several neurons, links and layers; 
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1 FIGs. 16A-C are graphs illustrating how different smoothing parameters are 

2 used in smoothing a trained GRNN; 

3 FIG. 1 7 shows the neural network architecture of the GRNN algorithm of the 

4 present invention; 

5 FIG. 1 8 shows the sensitivity of UA when used to detect flow rate reductions; 

6 FIG. 19 shows the sensitivity of the entering/leaving water temperature 

7 difference when used to detect flow rate reductions; 

8 FIG. 20 shows the sensitivity of UA when used to fouling of the heat 

9 exchanger tubes; 

1 0 FIG. 2 1 shows the sensitivity of APPR when used to detect fouling of the heat 

1 1 exchanger tubes; 

12 FIG. 22 shows the sensitivity of 7J isen(ropic when used to detect compressor faults; 

13 FIG. 23 shows the sensitivity of rj motor when detecting motor or transmission 

14 faults. 

1 5 FIG. 24 shows the sensitivity of UA C , and APPR C when detecting fouling faults; 

16 FIG. 25 shows the sensitivity of compressor isentropic efficiency when 

1 7 detecting compressor faults; 

1 8 FIG. 26 shows the sensitivity of motor efficiency when detecting motor faults; 

19 and 

20 FIG. 27 is a chart summarizing the relationships illustrated in FIGs. 18-26. 
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1 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

2 The inventor of the present invention discovered that it is difficult to detect and 

3 diagnose a fault in a chiller subsystem from raw data collected from sensors. Notably, it is 

4 difficult to set a meaningful fault threshold which is applicable to different chill subsystem 

5 implementations, i.e., is not customized for a particular implementation. To overcome this 

6 problem, the inventor of the present invention utilize a series of characteristic quantities 

7 which are calculated from the raw data using well known thermodynamic principles. The 

8 use of such characteristic quantities eliminates the need to individually tune the system of the 

9 present invention for any given HVAC installation. 

10 Importantly, the inventor of the present invention discovered that some 

1 1 characteristic quantities (CQs) are more sensitive than others in detecting certain faults. 

12 Table I below, summarizes the outputs (characteristic quantities) and the inputs used to 

13 determine them. 

14 The sensitivity of each CQ value to a given fault was determined by simulating 

1 5 various faults in a controlled experiment. Each fault condition in addition to a base case (no- 

1 6 fault) condition was simulated. Table II below pairs particular fault conditions with potential 

17 CQs that may be used for their detection. The inventor of the present invention developed 

1 8 a system for detecting and diagnosing faults utilizing the relationship summarized in Table 

19 II. 



10 



CA 02344908 2001-04-23 



TABLE I 





(Raw Data) 


Outputs (Characteristic Quantities) 


1 


T 

chws 


chilled water supply 
temp. (°F) 


APPR C 


evaporator heat exchanger approach 
(°F) 


2 


T 

chwr 


chilled water return 
temp. (°F) 


APPR C 


condenser heat exchanger approach 
(°F) 


3 


T 

cws 


condenser water 
supply temp. (°F) 


CHWTD 


chilled water temp, difference (°F) 


4 


m chw 


chilled water flow rate 
(gpm) 


CWTD 


condenser water temp, difference (°F) 


5 


m cw 


condenser water flow 
rate (gpm) 


UA C 


condenser conductance-area product 
(Btu/hr-°F) 


6 


T 

evap 


evaporator saturation 
temp. (°F) 


UAe 


evaporator conductance-area product 
(Btu/hr-°F) 


7 


T 

A cond 


condenser saturation 
temp. (°F) 


'Hisentropic 


isentropic heat transfer coefficient 
(Btu/hr-ft 2 - 0 F) 


8 


T 2 


compressor discharge 
temp. (°F) 


*H motor 


motor/transmission heat transfer 
coefficient (Btu/hr-fi 2 -°F) 


9 


P 


compressor power 


COP ovcraU 


overall coefficient of performance 


10 






COP 


compressor coefficient of performance 



12 


TABLE U 




13 


Fault Type 


Sensitive CQs 


14 


(1) Condenser water flow reduction 


UA C , APPR C , CWTD 


15 


(2) Chilled water flow reduction 


UA C , APPR C , CHWTD 


16 


(3) Evaporator tube fouling 


UA e and APPR e 


17 


(4) Condenser tube fouling 


UA C and APPR C 


18 


(5) Compressor internal faults 


P, COP como and rj^^ 


19 


(6) Motor/Transmission faults 


P, COP overaH and rj motor 



20 FIG. 1 is a block diagram of the overall structure of the model-based fault 

2 1 detection and diagnosis methodology of the present invention. The fundamental components 
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1 along with their links to other components and flow of information for fault detection and 

2 diagnosis are depicted, 

3 There are two sources of data used in this methodology, operating CQ's 

v 

4 calculated on-the-fly, and base-case CQ's which are interpolated from data stored in a look- 

5 up table. As shown in FIG. 1, measured inputs (and forcing functions), such as the chilled 

6 water supply temperature of the chiller and cooler water are fed to both the thermodynamic 

7 preprocessor 100 and the interpolator 102. The thermodynamic preprocessor 100 uses the 

8 measured inputs to generate the characteristic quantities (CQs). The operating CQ's are 

9 compared with base-case CQ's which are interpolated from empirical data stored in a base 

10 case lookup table 104. 

11 The look-up table 104 stores values of CQ's that were generated during a 

12 period in which fault-free operation of the chiller is assumed. The interpolator 102 

13 interpolates normal (base-case) CQ's from the data stored in the look-up table 104. The 

14 process for obtaining the operating and base case CQ values will be explained in detail 

15 below. 

16 The streams of operating and base-case CQ data are compared and errors 

1 7 (residuals) are computed. The residuals are fed to the Fault Classifier section 1 08 where the 

18 appropriate statistical analyses are employed to determine the significance of these errors. 

19 If they are indeed statistically significant, a fault detection alarm is produced. Otherwise, 

20 monitoring continues with the assumption that the chiller process is operating normally and 

21 that no faults are present. 

12 
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If faults are detected, the diagnostic classifier determines the type of fault and 
its likely location within the chiller. The CQs used in the methodology of the present 
invention were specifically selected to make isolating the fault location easier. As noted 
above, the inventor of the present invention discovered that certain CQs are sensitive to some 
faults and not to others. Thus, the careful selection of CQs makes it possible to distinguish 
between fault types. 

Once faults are detected and diagnosed, they are subsequently evaluated to 
determine the most suitable action to take. Fault Evaluation deals with comfort, safety, 
environment and economic issues. Action could be taken whenever: (1) comfort cannot be 
maintained, (2) equipment or personal safety is compromised, (3) environmental damage 
occurs (e.g. refrigerant leakage), or (4) reduced energy costs justify the service expense. 

In order to conduct on-line FDD using day-to-day operating data, it is desirable 
to have a model that can generate CQ's without the need for any manual intervention. The 
thermodynamic preprocessor 100 requires only eight measurement variables to generate the 
necessary CQ's immediately. See FIG. 2. Sets of CQ's are generated at two different time 
periods. As will be explained below, plural first sets of CQ variables (base-case CQ's), 
representative of normal (fault- free) system operating conditions, are generated and stored 
in the look-up table 104 during a preliminary period when a sufficiently wide range of loads 
is covered. A second set of CQ data (operating CQ f s) is generated "on-the-fly", i.e., 
generated contemporaneously with the receipt of the measured inputs to reflect the existing 
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1 state of the system, after the preliminary period has expired, and is compared with the first 

2 set of CQ variables to determine if faults have occurred. 

3 Both sets of CQ variables are generated in the same manner, the only 

4 difference being the time period in which each set is generated. Namely, the plural sets of 

5 base-case CQ's are generated during a preliminary period in which fault-free operation of the 

6 chiller is carefully controlled. The operating set of CQ f s are generated on-the-fly after the 

7 base-case CQ sets have been stored in the look-up table 104. 

8 The thermodynamic preprocessor 100 generates the operating CQ's using a 

9 thermodynamic data reduction program which was written in EES, an Engineering Equation 

10 Solver Software developed by S.A. Klein and FX. Alvarado, F-Chart Software, Inc. The 

1 1 EES program is particularly suited for solving thermodynamic equations; however, one of 

12 ordinary skill in the art will appreciate that other programming languages or software 

13 packages may also be used. 

1 4 The thermodynamic data reduction program consists of two procedures (FIGs. 

15 3 and 4) and equation sets for the condenser side input data (FIG. 5), evaporator side input 

16 data (FIG. 6), compressor input data (FIG. 7), evaporator equations (FIG. 8), compressor 

17 equations (FIG. 9), condenser equations (FIG. 10), expansion valve equations (FIG. 11) and 

18 for the thermodynamic properties and conversion constants (FIG. 12). It should be noted that 

19 the thermodynamic data reduction program uses general thermodynamic principles, and is 

20 not specifically tuned for a specific chiller subsystem. For this reason, the CQ f s have general 

21 applicability to all chiller subsystems. 

14 
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1 A procedure statement heads each procedure (FIGs. 3 and 4), and identifies the 

2 name and arguments of the procedure. Thus, the statement "Procedure LMTDj(T_chws, 

3 T_chwr, T_evap: LMTDe, T_we, ln_arg_e) ,f defines a procedure LMTD_e and a series of 

4 arguments (inputs and outputs). 

5 Table I, above, lists the input/outputs of the thermodynamic data reduction 

6 program. Moreover, FIG. 13 is a schematic of a chiller showing the location of each of the 

7 sensors used to obtain the measured inputs. 

8 The procedures (FIGs. 3 and 4) utilize function calls in which the function 

9 max(x,y) returns the larger of the two arguments (x, y). Thus, max(0, -2) returns "0". 

10 Correspondingly, the function ln(x) returns the natural log of argument x. 

1 1 The procedure LMTD_e (FIG. 3) is used to determine the logarithmic mean 

12 temperature difference between the evaporator saturation temperature (Tcvap), the chilled 

1 3 water supply temperature (T_chws), and the chiller water return temperature (T chwr) . The 

14 built-in EES error procedure halts the calculations if T_evap is greater than T_chws or if 

1 5 T chws is greater than T_chwr. 

16 The procedure, LMTD_c (FIG. 4) is used to determine the logarithmic mean 

17 temperature difference between the condenser saturation temperature (T_cond), the 

18 condenser water supply temperature (T_cws) and the condenser water return temperature 

19 (T_cwr). The built-in EES error procedure halts the calculations if T_cond is less than 

20 T cwr or if T cwr is less than T cws. 
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1 In developing the base case look-up table, measurement data for a wide variety 

2 of base-case operating conditions is compiled in a table which is accessed by the 

3 thermodynamic data reduction program. As will be appreciated by one of ordinary skill in 

4 the art, it is advantageous to collect measurement data for as wide a variety of operating 

5 conditions as possible in order to improve the accuracy of the interpolated CQ values. For 

6 example, it is desirable to calculate base-case CQ values for various times of day for both 

7 workdays and non- workdays, different weather conditions (ambient temperature, sunlight, 

8 winds, humidity), etc to fully reflect the wide range of heating/cooling loads encountered by 

9 the system. 

10 In FIGs. 5-7, the statement lookup( x, y) "F M is used to access entry "x" in the 

11 lookup table and assign the returned value to variable "y". The units for each variable 

12 accessed from the lookup tabic is provided in the comments following each lookup 

13 statement. Thus, in the above example, the units are in degrees Fahrenheit. 

1 4 As noted above, a first set of CQ variables, representative of normal (fault- free) 

15 system operating conditions, is generated during a preliminary period when a sufficiently 

1 6 wide range of loads is covered. It is assumed that sufficient data are collected and that there 

17 are no major unforeseen systematic changes during data collection. 

18 Fault detection is performed by comparing each respective operating CQ 

1 9 variable with a corresponding base-case CQ variable. Each set of CQ variables stored in the 

20 lookup table 104 includes at least eight (of the nine) input and ten (of the output variables 

21 listed in Table L 

16 
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As will be appreciated by one of ordinary skill in the art, there are many 
different combinations of input variables, whereas the look-up table only contains a relatively 
small subset of the possible combinations. Consequently, it is unlikely that a given set of 
eight measured inputs acquired on-the-fly will directly match one of the sets of input 
variables stored in the look-up table 104. Accordingly, the interpolator 102 determines a set 
of base-case CQ variables through a process of interpolation. 

The inventor of the present invention investigated a number of different 
interpolation schemes. Linear interpolation and regression against power form relations were 
found to be inaccurate and not general. According to a preferred embodiment, Neural 
Network architecture was found to be well suited to the interpolation. However, one of 
ordinary skill in the art will appreciate that other interpolation methods may be substituted. 

Neural network models are inspired by the human thought processes. Over 1 00 
billion biological neurons are present in the human brain. The connections between these 
neurons are called synapses and when the brain learns their strength is modified. 

Analogously, artificial neural networks contain artificial neurons that are 
connected via one-way information conduits called links. Weights are associated with these 
links that control the magnitude of the input signal entering the artificial neuron. These link 
weights simulate the physical and neuro-chemical characteristics of the biological synapse. 
Each artificial neuron in the network functions by first summing its scaled inputs and then 
applying a non-linear function to this sum to generate an output. FIG. 14 illustrates the 
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1 simple structure of a single neuron unit and FIG. 15 shows the general architecture of a 

2 interconnected network with several neurons, links and layers. 

3 A layer is a grouping of neurons and are of three basic types. The input layer 

v 

4 receives data from outside of the network, the output layer holds the network's final 

5 computational results and any layers in between these two are called hidden layers. 

6 A basic feature of neural networks is that they learn by example and are never 

7 programmed. Therefore, in order to train a network, it must be presented with a data sample 

8 that contains the inputs and their corresponding outputs. By an iterative process, the network 

9 gradually learns the input-output relationship and then, depending on the application, is used 

10 for prediction, correction or pattern recognition. Once trained, they can resolve numeric 

1 1 problems that are otherwise solved by conventional regression analysis. In the network, the 

1 2 inputs are equivalent to the independent variables used in regression, the dependent variables 

13 become the outputs and the observations are the sample patterns used for training the 

14 network. 

1 5 Neural networks fall into the two basic categories of: supervised networks and 

16 unsupervised networks. Supervised networks are those that make predictions, classify 

17 patterns or make decisions based on a set of learned inputs and outputs. Some examples of 

1 8 these are the Probabilistic Neural Networks (PNN), the General Regression Neural Networks 

19 (GRNN) and Back Propagation Networks. These kinds of networks are trained to make 

20 predictions, classifications, or decisions when given a large enough number of accurate 

21 classifications or predictions from which it can leam. 

18 



CA 02344908 2001-04-23 



1 Unsupervised networks, without being taught how to categorize, are able to 

2 classify a given number of input and output patterns into a specified number of categories. 

3 They do this by clustering the training patterns using their proximity in n dimensional space 

4 where n is the number of inputs. The network usually clusters the data into the maximum 

5 number of categories it is presented with. Kohonen Networks are examples of this network 

6 type. 

7 The GRNN is one type of neural network that is well suited to interpolation. 

8 It is based on the estimation of a probability density function of a vector random variable, 

9 X, and a scalar random variable, Y. If the joint probability density function of these variables 

1 0 are both known then the conditional probability density function and the expected value can 

1 1 be computed. The estimated value of Y for a given X is presented in the following general 

12 regression equation: 

OO 

\Yf(X,Y)dy 

13 E[Y\X) = ^ (1) 

\f{X,Y)dy 

14 where 

15 J^Y\ X] = conditional mean of Y on X 

1 6 f(X* Y) = known joint continuous probability density function 
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The probability density function JpC Y), is estimated from sample observations 
of X and Y when it is unknown. For a non-parametric estimate of f(X,Y), the Parzen 

estimation f\X 9 Y), is used by the GRNN. It is defined by the following equation for the 
observed sample observations, X t and Y f of the vector X and scalar Y: 



where 



and 



1 



1 



f'(X,Y) = - x -Zfxfy 



(2) 



fx = ex P 



(x-x^jx-x,) 

2a 2 



(3) 



/ K =exp 



2a 2 



(4) 



and 



= the number of sample observations 
= the dimension of the X vector 
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cr — the standard deviation (or smoothing parameter) 

An estimate for the desired mean of Fat any given Xis derived in Equation (5) 
by combining Equations (1) and (2) and performing the integration after first interchanging 
the integration and summation operation. 



Y(X) = 



( D 2 ^ 
K 2a J 



£exp 



' D 2 ^ 



1 = 1 



v 2 ^y 



(5) 



where, the scalar function D i is given by 



=(x-x,Y{x-x,) 



(6) 



The main algorithm of the GRNN model is expressed by Equations (5) and (6). 

The estimate Y(X) is a weighted average of all the observed samples, Y i9 where each sample 

is weighted in an exponential manner according to the Euclidean distance, D g9 from each^. 
This appropriate weighting is explained by the inversely proportional relationship between 



the expression exp 



andD,. That is, as jD ; increases, exp 



' A 2 



2a 2 . 



decreases and vice- 



versa. An optimum value for the smoothing parameter, cr, is determined using iterative or 
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genetic adaptive methods. Iterative methods such as the Holdout method and the Wiggle 
method are well known by those of ordinary skill in the art, and do not form part of the 
present invention. Larger values of cr improve smoothness of the regression surface. 



Notably, cr must be greater than 0 and can usually range from .01 to 1 with good results. 

FIGs. 16A-C illustrates how different smoothing parameters are used in 
smoothing a trained GRNN. In this example, one input value is used to predict one output 
value. As shown in the graph, of the forty input and output patterns only input patterns 10, 
20, and 30 produce an output value of 1. 



Equations (5) and (6). The Euclidean distance, £>,., is computed by the links between the 
input layer and the first hidden layer. Based on observed samples, X i9 and smoothing 
parameter, a y the expression, 



is computed. A node in the second hidden layer takes the sum of the exponential values of 
all samples. In other nodes of this same layer, the products of the exponential values and the 
corresponding observed Y t for each sample observation are computed. The node in the third 
hidden layer computes the sum of all these product values, which is then supplied to the 
output node where the ratio between it and the previous sum is calculated. 



FIG. 17 represents the neural network architecture of the GRNN algorithm of 
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1 The suitability of the GRNN, however, is attributed to several important 

2 features and makes it convenient for online implementation. In the GRNN network only a 

3 single parameter is estimated. Unlike other networks, a once through, non-iterative training 

4 process with a highly parallel structure is involved. Compared to conventional regression 

5 analysis, the specification of the underlying regression function, bounds of the independent 

6 variables, initial convergence values and convergence criteria are not required beforehand. 

7 Additionally, the algorithm provides smooth transitions from one observed value to another 

8 even with sparse and noisy data in a multidimensional measurement space and can be used 

9 for any regression problem where an assumption of linearity is not justified. 

10 Since the thermodynamic and the GRNN models involve functions and 

1 1 variables in a multi-dimensional space, it is important to know what the truly independent 

12 measurement variables are. Here, the term independent refers to those variables that are 

13 forcing functions of the chiller process. Information about variable independence is 

14 necessary for understanding the available data, achieving greater insight to data 

1 5 manipulation, accurately assigning inputs and outputs, and successfully training the GRNN 

16 model to be able to predict normal system behavior. For the chiller process, only five 

17 measurement variables may be considered independent. This number is determined from 

1 8 studying the various equations that define the chiller process more closely in the discussion 

19 that follows. 

20 An overall energy balance of the entire refrigeration system is: 
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Power — Q ~Q . 

2 Ocond T ^cvf^pJ^cwr ^cms) ^g^ 

3 £L/ =^V(^2 -Ki^^^cand (9) 

4 The enthalpy difference, Ah ref cond , between the entering and leaving refrigerant 

5 in the condenser is a function of the condenser saturation temperature, T cond : 

g ^'nf, ootid ~ f^cotid ) (10) 

8 ^=™^K~\) = <f^^ (12) 

9 Here, the enthalpy difference, Ari retevQpp between the entering and leaving 
10 refrigerant in the evaporator is a function of the evaporator saturation temperature, T evap : 

=X T e« P ) (13) 

12 In addition to the above relations, the following relations for the chiller Power 

13 and Capacity, are available from chiller manufacturers' experiments: 
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, Powsr=f(T emp ,T caxi ) (l4) 

2 Qevap f^^extp^cond} (15) 

3 Equations (7) to (15) fully define the thermodynamic data reduction model. 

4 An inventory of the number of equations and the number of unknowns yield nine and 

5 fourteen, respectively. Thus the difference between these two numbers tells us that five 

6 variables must be independently measured for these model equations to work. 

7 The inventor of the present invention selected the following five independent 

8 measurement points condenser supply temperature, T cws , the chilled water supply 

9 temperature, T chws , the chilled water return temperature, T chwr9 the condenser water flow rate, 

10 ff\ W9 and the chilled water flow rate, However, the selection of other independent 

1 1 measurement points is contemplated and falls within the scope of the present invention. 

12 The condenser water supply temperature, T cws and the chilled water return 

13 temperature, T chwr are the entering water temperatures to the chiller subsystem. The chiller 

14 subsystem has no direct effect on the magnitudes of these two water streams until after they 

1 5 cross the system boundary. Therefore, the temperatures, T chwr and T cws are independent. They 

16 represent only a part of the building load and heat rejection, respectively. The flow rates, 

17 *hhw an d JVnv, are additional required independent variables. The final measurement 
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1 variable, the chilled water supply temperature, T chwsy is termed independent because its 

2 magnitude is associated with a control set-point. 

3 The essence of the Fault Detection scheme relies on the determination that a 

4 fault truly exists based on the deviation of a Characteristic Quantity (CQ) exceeding some 

5 predetermined threshold or criterion. This deviation is measured from an established base 

6 case representing normal operating conditions. 

7 As described above, the CQ values are computed from measurement data 

8 gathered from sensors. Due to the many perturbations that exist, there is an inevitable chain 

9 of causes and effects that occur and form the basis for uncertainties in the measurement data. 

10 Some of the factors that may cause perturbations are power supply variance, hysteresis, 

1 1 precision and bias errors and general instability. These factors contribute to the accuracy 

12 ratings of the measurement sensors and therefore, are factors that contribute to the 

13 uncertainties in CQ deviations. These uncertainties in CQs determine the threshold and thus 

14 the sensitivity for fault detection. A tighter uncertainty region corresponds to more sensitive 

15 detection. 

16 Numerous situations involve the computation of an important quantity, 7, 

17 instead of its direct measurement. If such a quantity is determined from Another directly 

1 8 measured quantities X Jy X 2 , . . . , X N , a functional relation between these quantities and 7 may 

19 be expressed as Y z=z f(X x ,X 2 ,„.,X N ). This function f not only expresses a physical law but 
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a measurement process. In fact, all the quantities that could contribute a significant 
uncertainty to the measurement results are accounted for by / 



Furthermore, if ^represents the estimate of the output quantity, 7, and x l5 x 2 , 
. . . , x N are input estimates for the values of the N input measurements, X j9 X>, . . . , X N9 then Y 

can be expressed as Y=f(pc l9 x l9 ... 9 x N ). In order to obtain the combined standard 
uncertainty, u c (Y) of the computation result Y (an estimated standard deviation of the 

result), a first-order Taylor series approximation of Y=f(X X9 X 29 .„,X N ) is performed. This 
gives the following expression known as the Law of Propagation of Uncertainty: 



N f jr\ AM N <or jqr 

(=1 \6&iJ r =l yW+l^i 



-|]/2 



(16) 



where 



= sensitivity coefficients 



Xj~Xf 
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t^Xi) = standard uncertainty associated with the input estimate, x { 

n(x i 9 Xj ) == estimated covariance associated with x i and Xj 

If the individual measurement inputs are assumed uncorrected and 
random, then the covariance term is zero and U C (Y) is simplified to: 



u c (Y) = 



2 nl/2 



(17) 



By manner of example, two CQs that are sensitive to tube fouling in the 
condenser arc its approach, APPR ci and conductance-area product, UA C . Using the 
simplified form of the uncertainty propagation law in Equation (17), expressions can 
be derived for the uncertainties in UA C and APPR C . 

Equation (18) that defines the conductance-area product of the 
condenser, UA C , in terms of four independently measured variables and a fluid 
property, C pw . This property is the specific heat of water and actually varies with 
temperature conditions but was assumed constant for this application since the range 
of variation is insignificant. 



UA C = >?w C 



pw 



T —T ^ 

COtvi CMS 



(18) 
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The UA C is a f(Tcw^cw^cond> J Kw)* Equation (16) was used to express the 



most probable uncertainty in UA C as: 



u c {UA c ) = 



#JA C 



(sua. 



J 



f SUA, 



V^Lr ) \STcond ) 



SUA C 



1/2 



4 

5 



where, 



(19) 



(20) 



SUA, 



mC. 



^cwr i^cond ^cwr) 



(21) 



Sh„.. 



(^cond ^cws 
(Tcond ~ ^vir)> 



(22) 



(23) 



10 



and, i(T^\l(T^\UiT a ^) and are the uncertainties of the measured 



11 
12 
13 



condenser water supply temperature, condenser water return temperature, refrigerant 
condenser saturation temperature and the water mass flow rate, respectively. These 
uncertainties are equivalent to the accuracy of the measurement sensors. 
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The approach, APPR C is a f(T cwr , T con ^ . The most probable uncertainty 



in APPR C may be expressed as: 



u{APPR c )-. 



mPPR c 



mppR„ 



cond 



1/2 



(24) 



where 



5 
6 



and 



= -1 



(25) 



<3iPPR c 



+1 



(26) 
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For the threshold analysis, only precision errors were used. The 
following sensor precision errors, developed using time series analysis, were used to 
evaluate the effect of error on CQ: 

•Temperature sensor: ±0.3°F and ±0.6°F 

•Water flow rate sensor: ±5% 



•Power kW-meter: 



±10% and ±20% 
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1 Fault detection thresholds are illustrated in FIGs. 18 -24. Two types 

2 of horizontal lines are used in defining these thresholds. The no-fault line represents 

3 no flow rate reduction (or increase), no fouling and no inefficiencies in either the 

4 compressor or the motor. The critical line was determined by the following 3 -step 

5 procedure: 

6 1 . Locate the points where the uncertainty boundaries intersect the no-fault 

7 line. 

8 2. Project these points vertically to the points that intersect the 0.0°F curve 

9 (shown as black dots; for ±0.3 °F, refer to FIG. 1 8). 

10 3 Draw horizontal lines, parallel to the abscissa, through these intersection 

11 points. 

12 These horizontal lines are the critical lines and they define a region 

1 3 (shown as hatched in FIG. 1 8) within which it is possible for faults to go undetected. 

14 Beyond these critical lines (or thresholds), however, faults may be detected with more 

15 confidence despite the errors present in the measurement sensors. 

16 For example, water flow rate abnormalities (reduction or increase) are 

17 detected most effectively by deviations in the conductance-area product and the 

1 8 entering/leaving water temperature difference of the given heat exchanger (evaporator 

19 or condenser). FIG. 18 shows the sensitivity of UA when used to detect flow rate 

20 reductions. The heavy solid curve represents the UA deviation in the face of no errors 

21 in the temperature and flow sensors (i.e., 0.0°F/0% of full scale, gpm). The dashed 
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1 and light solid curves represent the UA deviation due to errors of ±0.3°F/±5% gpm 

2 and ±0.6°F/±5% gpm, respectively. 

3 Here, the critical lines are drawn using the uncertainty boundaries 

4 defined by ±0.6°F/±5% gpm. Based on a ±0.6°¥/±5% gpm precision of these sensors, 

5 it can be concluded that UA deviations between —23.6% and +23.6% would not 

6 indicate a fault. Flow rate reductions greater than 41.0% and flow rate increases 

7 greater than 68.0% would be needed to indicate faults. 

8 FIG. 1 9 is similar to FIG. 1 8 in terms of the precision errors represented 

9 by the heavy, dashed and light curves. However, the characteristic quantity used here 

10 is the entering/leaving water temperature difference instead of the conductance-area 

1 1 product. It is apparent from the critical lines drawn in FIG. 1 9 that the thresholds for 

1 2 detecting flow rate reduction is lower than those depicted in FIG. 1 8. In others words, 

13 this characteristic quantity has a greater fault detection sensitivity. 

14 Quantitatively, here, it is concluded that water temperature difference 

15 deviations beyond -14.3% and +14.3% are used to detect flow rate increases greater 

16 than 16.0% and flow rate reductions greater than 13.0%, respectively. 

17 Fouling describes a decrease in the effectiveness of the heat exchanger 

1 8 tubes, and is detected most effectively by deviations in the conductance-area product 

1 9 and the approach of the given device. FIG. 20 shows the sensitivity of UA when used 

20 to detect fouling. A fouling factor of 1 represents no fault and in this study, is the 

21 normalized equivalent of 0.00075 hr-ft 2 /BTU. A factor greater than 1 is an increase 
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1 in fouling and lesser than 1 is a decrease. The critical lines indicate that UA 

2 deviations beyond -19.6% and +19.6% are used to detect fouling factors greater than 

3 1.6 and lesser than 0.3, respectively. 

v 

4 FIG. 21 shows the sensitivity of the Approach when used to detect 

5 fouling. The critical lines indicates thatAPPR deviations beyond-1 7.6% and +1 7.6% 

6 are used to detect fouling factors lesser than 0.7 and greater than 1.3, respectively. 

7 These numbers, compared to those for UA, suggest that the APPR has lower fouling 

8 detection thresholds. 

9 Some internal faults of the compressor are detected most effectively by 

10 deviations in its isentropic efficiency, h iseniropic . FIG. 22 shows the sensitivity of 

1 1 Visentropic when used to detect faults. A compressor fault factor of one represents no 

12 fault. In contrast to the fouling factor, a compressor fault factor greater than one is 

13 an improvement of the fault condition and a factor lesser than one is a worsening of 

14 the situation. The critical lines indicate that rj isentropic deviations beyond -2.1% and 

15 +2.1% are used to detect fault factors lesser than 0.8 and greater than 1 .2 respectively. 

1 6 Deviations in the motor efficiency, rj motot > are used to detect faults in the 

1 7 motor or transmission housing. FIG. 23 shows the sensitivity of rj motor when detecting 

18 motor or transmission faults. The heavy solid curve represents the 7] mo(or deviation 

19 when the kW-meter is error free (i.e., 0% of full scale, kW). The dashed and light 

20 solid curves represent the rj motor deviation due to errors of ±10% kW and ±20% kW, 
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1 respectively- A motor fault factor of one represents no fault and in this study, is 

2 equivalent to 0. 85 . However, factors greater than one represent worsening conditions . 

3 Here, the critical lines are drawn using the uncertainty boundaries 

4 defined by the ±20% kW error. Based on this kW-meter precision, it is concluded 

5 that rj motor deviations less than -20% is used to detect fault factors greater than 1 .3. 

6 Equations (27)-(29) below are similar in form to the previous three 

7 relations for the condenser, and are used to model the chilled water reduction. Here, 

8 the chilled water flow rate reduction causes an increase in the return temperature and 

9 a decrease in the evaporator temperature for a given load and chilled water set point 

10 temperature. 

11 'K^fouit = m^Jl - (27) 

1 2 Q„ w - C jnv (T dM ,. - T dn ^ ) (28) 

h «= U7 { NETUBE j < 29 > 

1 4 where, ^dxwjhiM , factor me , NEPASS, NETUBE and D e are the chilled water flow due to 

15 the flow reduction fault, the percentage reduction in chilled water flow rate, the 

16 number of evaporator tube passes, the number of evaporator tubes and the diameter 

17 of the evaporator tubes, respectively. 
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1 The CQs UA €9 APPR e and CHWTD show the most significant deviations 

2 from the no-fault base-case. The CQ relationships here are similar to those used for 

3 the condenser water reduction fault. CHWTD is the most sensitive and its sensitivity 

4 is unaffected by load changes. 

5 Equations (30) and (31) are used to model condenser tube fouling. 

6 There are some conflicting interpretations of the fouling factor in the literature. Here, 

7 fouling is defined as the resistance to heat transfer due to the tube material as well as 

8 to any thin film build-up on its surfaces. Equation (30) is used to simulate an increase 

9 (or decrease) of the level of fouling in the condenser tubes. It is apparent from 

10 Equation (31) that as fouling increases, UA C decreases. 

11 KjaUt^K- fc (30) 

12 UA C = j- (31) 

13 where, factor fc is the fault modification factor used to simulate degrees of fouling in 

1 4 condenser (unity = normal) and R c fi3uit is the resistance to heat transfer associated with 

15 the condenser tube material, including the fouling factor for faulty conditions. 

1 6 The UA C , and APPR C are CQs that show the more significant deviations 

17 from the no-fault base-case. FIG. 24 shows how these CQs deviate from the base- 

18 case condition as the fault condition changes (i.e., as fouling increases or decreases). 
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Similar to the flow reduction faults, the different CQs behave differently. For 
example, at a base-case load of 1500 tons, a 100% fouling increase by the condenser 
tubes corresponds to a 27% decrease and a 54% increase in UA C and APPR 
respectively. However, in contrast to condenser flow reduction faults, CWTD shows 
no significant change. This behavior is important in distinguishing the fouling fault 
from the flow reduction fault. 

At 500 tons, the same fouling condition corresponds to a 29% decrease 
and 63% increase in UA C and APPR C , respectively. It is apparent from these results 
that the more sensitive CQ is the APPR C and that varying loading conditions affect it 
more than they did UA C . 

Similar trends exist for a 50% fouling decrease condition. CQs that 
increased before decreased and vise- verse. Fouling decrease information is as useful 
as fouling increase information. Here, FDD may be used to check the integrity of 
scheduled maintenance, repair or cleaning of fouled tubes. 

Equations (32) and (33) are used to model evaporator tube fouling. 

Kjauit^K-fartWfe (32) 
it 



UA=- 



ei 



L. 1 ,p 

h rh fault 

"e,i r e"e,o (33) 
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1 where,factor fe , is the fault modification factor used to simulate degrees of fouling in 

2 the condenser (unity ■ normal) and R e and R eJault are the resistance to heat transfer 

3 associated with the evaporator tube material, including the fouling factor for normal 

4 and for faulty conditions, respectively. 

5 The UA e and APPR € are CQs that show the more significant deviations 

6 from the no-fault base-case. The CQ relationships here are similar to those used for 

7 the condenser tube fouling fault. APPR e is more sensitive and its sensitivity is slightly 

8 affected by load changes. 

9 The Weisner correlation, based on the measured performance of 

10 centrifugal compressors, is used for modeling compressor internal faults. The 

1 1 modification of the reference polytropic efficiency, rf refen in Equation (34) is used to 

12 simulate possible faulty conditions within the compressor. 

1 3 Vrefafault = Prefer f aCt O& mp (34) 

14 where, factor CQmp and rj refer/auln are the compressor fault modification factor, and the 

15 modified rj refer due to faulty conditions, respectively. 

1 6 The compressor isentropic efficiency, rj iseniropic is the CQ that shows the 

17 most significant deviations from the no-fault base-case. FIG. 25 illustrates a greater 

18 decrease in the isentropic efficiency as the compressor fault factor decreases. Based 

19 , on the manner in which this fault is simulated, a decrease in the fault factor means a 
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worsening compressor internal fault condition. Varying the loading conditions does 
not appear to have a great impact on the deviation of the isentropic efficiency. 

Whenever there is a fault in the motor or transmission system of the 
compressor, the overall power draw (P) is expected to increase compared to the actual 
power used for the compression process. The modification factor in Equation (35) 
simulates this increase. Equation (35) defines the combined efficiency of the motor 
and transmission. An energy balance on the compressor establishes the relation 
between the states of the refrigerant at the suction and discharge conditions. 

PfauU^P'fi^or^ (35) 

7 lno/or ~~ p (36) 
"fault 

where, P faul „ factor moior and A r/^^ are the compressor power consumption due to 
faulty conditions, the motor fault modification factor and the change in refrigerant 
enthalpy between compressor's suction and discharge states, respectively. 

The motor efficiency is the CQ that shows the most significant 
deviation from the no-fault base-case. FIG. 26 shows that there is a greater decrease 
in the motor efficiency as the motor fault factor increases (i.e., as the power draw 
increases for a given loading condition). From this simulation, the motor efficiency 
deviation (not the motor efficiency) is independent of loading conditions. That is, at 
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1 any loading condition, the percentage deviations due to an increased fault condition 

2 will show the same magnitudes. 

3 Refrigerant loss results from leakage in chiller subsystems, which has 

4 the net effect of reducing overall system performance. It is often confirmed by 

5 bubbles present in the refrigerant (as seen from the liquid line sight glass), frost 

6 formation, abnormal superheat and subcooling temperatures, and abnormal unit 

7 operating pressures. A temperature difference between the low pressure liquid line 

8 temperature (entering the evaporator) and the suction line temperature (entering the 

9 compressor) was computed and used as a characteristic quantity to flag the occurrence 

10 of refrigerant leaks . 

11 FIG. 27 summarizes the diagnostic patterns exhibited by the CQs 

12 discussed earlier. The symbols "0", and signify no significant change, an 

13 increase and a decrease compared to base-case values, respectively. In the Type 

14 column, "N" "F" and "N/A" stands for Normal, Fault and Not Applicable, 

15 respectively. Two normal operating conditions called "Load change** and "Heat- 

16 rejection change" are included. Respectively, these are quantities representing a 

17 change in the chiller load and a change in the amount of heat rejected by the 

1 8 condenser. They are modeled by varying the chilled water return and condenser water 

1 9 supply temperatures, respectively. These conditions are included to show that it is 

20 possible to distinguish them from fault conditions using the appropriate set of CQs. 
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Although, these CQs effectively isolate faults on the component level, 
the diagnostic patterns help to further distinguish between the various fault conditions 
that may exist in a particular component. For example, the increase and decrease in 
the condenser approach, APPR C and the condenser conductance-area product, UA C , 
respectively, would signify that there is a problem in the condenser and not in the 
other components of the chiller sub-system. However, this same CQ change pattern 
is both exhibited by an increase in condenser tube fouling and a decrease in the 
condenser water flow rate. To distinguish between these two conditions, FIG. 27 
shows that the condenser water flow rate decrease is accompanied by an increase in 
the condenser water temperature difference, CWTD, whereas the condenser tube 
fouling is not. 

For most performance monitoring studies, the COP overall is a CQ that is 
often used to determine significant performance changes. The COP overall is included 
in FIG. 27 to illustrate that although it gives some indication of changes within a 
chiller sub-system, by itself, it is quite vague in isolating faulty components or in 
diagnosing faults within a given component. Therefore, it must be used in 
combination with the other CQs, mentioned here, for useful diagnoses to be made. 

According to one modification of the above-described FDD 
methodology, the use of a base-case lookup table is eliminated. Notably, the use of 
parallel fault detection and diagnosis in which CQ residuals are derived from 
comparing a particular chiller to other identical (similar model, capacity and run-time) 
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chillers within the same HVAC system eliminates the need for a lookup table. In 
parallel FDD, the data from a designated chiller is considered to be the base-case, and 
faults are detected by comparing the data from another chiller with the data from the 
designated chiller. In so doing, faults in a particular chiller or chillers are determined 
over a predetermined time period. 

The parallel FDD methodology is substantially similar to the above 
described (serial) FDD methodology in all relevant respects. As will be appreciated 
by one of ordinary skill in the art, in parallel FDD the CQ's from the designated 
chiller and the other chiller(s) are determined on-the-fly, and the difference (residuals) 
between the two streams of CQ's is used to detect faults in the previously described 
manner. 

One of ordinary skill in the art will appreciate that both of the above- 
described embodiments of the present invention utilize CQ's and thus do not require 
modification or tuning of the thermodynamic model. Accordingly, the FDD system 
of the present invention may readily be included in an existing facility management 
system, provided that the appropriate sensors are already installed. Use of the first 
embodiment merely requires the collection and storing of base-case CQ's whereas 
base-case CQ's are calculated on the fly according to parallel FDD. 

Another benefit of the present invention relates to threshold used to 
detect the occurrence of faults. Both of the above-described embodiments detect of 
faults by comparing operating CQ's to base-case CQ's. A fault is detected when one 
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or more CQ's exceed a predetermined threshold which is determined in relation to the 
sensor accuracy. Accordingly, the false detection of faults is minimized. 

Yet another benefit of the claimed invention relates to the classification 
of faults. Notably, the inventor of the present invention discovered that certain CQ's 
are more sensitive than others to the occurrence of faults, and that faults may 
accurately be diagnosed by examination of which CQ's exceed the predetermined 
threshold. 

While various embodiments of the present invention have been shown 
and described, it should be understood that other modifications, substitutions and 
alternatives are apparent to one of ordinary skill in the art. Such modifications, 
substitutions and alternatives can be made without departing from the spirit and scope 
of the invention, which should be determined from the appended claims. 

Various features of the invention are set forth in the appended claims. 
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WHAT IS CLAIMED IS: 

1 1 . A fault detection system for an HVAC system including sensors for 

2 measuring the performance of one or more of a condenser, a compressor, an 

3 evaporator, and a chiller, said fault detection system comprising: 

4 a processing means for generating characteristic quantities (CQ's) from a 

5 plurality of measured inputs; 

6 means for storing plural CQ values generated by said processing means in an 

7 initial period in which fault-free operation of the HVAC system is assumed, CQ 

8 values being generated for each of plural different measured input values; 

9 means for producing interpolated base-case CQ values from said plural stored 
10 CQ values for a given set of measured inputs; 

1 ] means for detecting a fault when a difference between actual CQ values, 

1 2 calculated by said processing means using said given set of measured inputs, and said 

13 interpolated base-case CQ values varies from a predetermined threshold value; 

14 means for classifying a detected fault based on which ones of said actual CQ 

15 values varies from said interpolated CQ values. 

1 2. The fault detection system according to claim 1 wherein said 

2 processing means uses measured inputs selected from the group comprising: chilled 

3 water supply temperature; chilled water return temperature; condenser water supply 

4 temperature; chilled water flow rate; condenser water flow rate; evaporator saturation 
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5 temperature; condenser saturation temperature; and compressor discharge 

6 temperature. 

1 3. The fault detection system according to claim 1 wherein said 

2 interpolator is a General Regression Neural Network. 

1 4. The fault detection system according to claim 1 wherein said 

2 interpolator is a Probabilistic Neural Network. 

1 5. The fault detection system according to claim 1 wherein said 

2 interpolator is a Back Propagation Network. 

1 6. The fault detection system according to claim 1 wherein said CQ's are 

2 selected from the group comprising: evaporator heat exchanger approach; condenser 

3 heat exchanger approach; chilled water temperature difference; condenser water 

4 temperature; condenser conductance-area product; condenser conductance-area 

5 product; isentropic efficiency; motor/transmission efficiency; overall coefficient of 

6 performance; and compressor coefficient of performance. 

1 7. The fault detection system according to claim 6 wherein condenser 

2 water flow rate abnormalities are diagnosed when one of said condenser conductance- 
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3 area product, said condenser heat exchanger approach, and said condenser water 

4 temperature difference exceeds said predetermined threshold value for said CQ. 

1 8. The fault detection system according to claim 6 wherein chilled water 

2 flow rate abnormalities are diagnosed when one of said evaporator conductance-area 

3 product, said evaporator heat exchanger approach, and said chilled water temperature 

4 difference exceeds said predetermined threshold value for said CQ. 

1 9. The fault detection system according to claim 6 wherein chilled water 

2 flow rate abnormalities are diagnosed when one of said evaporator conductance-area 

3 product, said evaporator heat exchanger approach, and said chilled water temperature 

4 difference exceeds said predetermined threshold value for said CQ. 

1 10. The fault detection system according to claim 6 wherein evaporator 

2 tube fouling is diagnosed when one of said evaporator conductance-area product and 

3 said evaporator heat exchanger approach exceeds said predetermined threshold value 

4 for said CQ. 

1 11. The fault detection system according to claim 6 wherein condenser tube 

2 fouling is diagnosed when one of said condenser conductance-area product and said 
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condenser heat exchanger approach exceeds said predetermined threshold value for 
said CQ. 

12. The fault detection system according to claim 6 wherein an internal 
compressor internal fault is diagnosed when one of said compressor coefficient of 
performance, said isentropic efficiency, and compressor power draw exceeds said 
predetermined threshold value for said CQ. 

13. The fault detection system according to claim 6 wherein a motor 
transmission fault is diagnosed when one of said overall coefficient of performance, 
said motor/transmission efficiency, and motor power draw exceeds said 
predetermined threshold value for said CQ. 

14. The fault detection system according to claim 6 wherein said 
predetermined threshold includes upper and lower critical value levels for each said 
characteristic quantity, said upper and lower critical value levels being determined in 
accordance with a rated precision error of a sensor used to acquire a selected said 
measured input used to calculate said characteristic quantities. 

15. A method for detecting faults in a chiller subsystem of a facility cooling 
system, comprising: 
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3 providing plural base data sets of measured inputs, each said set of measured 

4 inputs including sensor data for a different base operating condition of the chiller; 

5 computing a set of characteristic quantities for each of said plural base data 

6 sets; 

7 storing said plural sets of characteristic quantities in a memory means; 

8 providing a test data set for a test operating condition of the chiller; 

9 computing a set of characteristic quantities from said test data set; 

1 0 producing an interpolated base set of characteristic quantities from among said 

1 1 plural sets of characteristic quantities stored in said memory means using said test 

12 data set; 

13 detecting a fault when at least one characteristic quantity exceeds a 

1 4 predetermined threshold range; 

15 diagnosing a fault in relation to which ones of said characteristic quantities 

1 6 exceed said predetermined threshold range. 

1 16. The method according to claim 15 wherein said predetermined 

2 threshold range includes upper and lower critical value levels for each said 

3 characteristic quantity, said upper and lower critical value levels being determined in 

4 accordance with a rated precision error of a sensor used to acquire data items 

5 contained in said test data set. 
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17. The method according to claim 15 wherein said measured inputs are 
selected from the group comprising: {chilled water supply temperature; chilled water 
return temperature; condenser water supply temperature; chilled water flow rate; 
condenser water flow rate; evaporator saturation temperature; condenser saturation 
temperature; and compressor discharge temperature. 

1 8 . The method according to claim 1 5 wherein said characteristic quantities 
are selected from the group comprising: evaporator heat exchanger approach; 
condenser heat exchanger approach; chilled water temperature difference; condenser 
water temperature; condenser conductance-area product; condenser conductance-area 
product; isentropic efficiency; motor/transmission efficiency; overall coefficient of 
performance; and compressor coefficient of performance. 

19. The method according to claim 1 8 wherein condenser water flow rate 
abnormalities are diagnosed when one of said condenser conductance-area product, 
said condenser heat exchanger approach, and said condenser water temperature 
difference exceeds said predetermined threshold value for said CQ. 

20. The method according to claim 1 8 wherein chilled water flow rate 
abnormalities are diagnosed when one of said evaporator conductance-area product, 
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said evaporator heat exchanger approach, and said chilled water temperature 
difference exceeds said predetermined threshold value for said CQ. 

21. The method according to claim 18 wherein water flow rate 
abnormalities in said chilled water flow are detected when one of said evaporator 
conductance-area product, said evaporator heat exchanger approach, and said chilled 
water temperature difference exceeds said predetermined threshold value for said CQ. 

22. The method according to claim 18 wherein evaporator tube fouling is 
diagnosed when one of said evaporator conductance-area product and said evaporator 
heat exchanger approach exceeds said predetermined threshold value for said CQ. 

23. The method according to claim 18 wherein condenser tube fouling is 
diagnosed when one of said condenser conductance-area product and said condenser 
heat exchanger approach exceeds said predetermined threshold value for said CQ. 

24. The method according to claim 18 wherein an internal compressor 
internal fault is diagnosed when one of said compressor coefficient of performance, 
said isentropic efficiency, and compressor power draw exceeds said predetermined 
threshold value for said CQ. 
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25. The method according to claim 1 8 wherein a motor transmission fault 
is diagnosed when one of said overall coefficient of performance, said 
motor/transmission efficiency, and motor power draw exceeds said predetermined 
threshold value for said CQ. 

26. A fault detection and diagnosis system for a HVAC system including 
at least two chillers, each chiller being equipped with sensors for measuring the 
performance of one or more condenser, a compressor, an evaporator, and a chiller, 
said fault detection system comprising: 

processing means for calculating characteristic quantities (CQ's) from a 
plurality of measured inputs for each of the chillers, one chiller being designated as 
a base-case chiller; 

means for detecting a fault when a difference between CQ values for the base 
case chiller and CQ values for other ones of the chillers exceeds a predetermined 
threshold range for at least one of said CQ values; 

a fault classifier for classifying a detected fault based on which ones of said 
actual CQ values exceed said threshold range. 

27. The fault detection system according to claim 26 wherein said CQ's are 
selected from the group comprising: evaporator heat exchanger approach; condenser 
heat exchanger approach; chilled water temperature difference; condenser water 
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temperature; condenser conductance-area product; condenser conductance-area 
product; isentropic efficiency; motor/transmission efficiency; overall coefficient of 
performance; and compressor coefficient of performance. 

28. The fault detection system according to claim 27 wherein condenser 
water flow rate abnormalities are diagnosed when one of said condenser conductance- 
area product, said condenser heat exchanger approach, and said condenser water 
temperature difference exceeds said predetermined threshold value for said CQ. 

29. The fault detection system according to claim 27 wherein chilled water 
flow rate abnormalities are diagnosed when one of said evaporator conductance-area 
product, said evaporator heat exchanger approach, and said chilled water temperature 
difference exceeds said predetermined threshold value for said CQ. 

30. The fault detection system according to claim 27 wherein chilled water 
flow rate abnormalities are diagnosed when one of said evaporator conductance-area 
product, said evaporator heat exchanger approach, and said chilled water temperature 
difference exceeds said predetermined threshold value for said CQ. 

3 1 . The fault detection system according to claim 27 wherein evaporator 
tube fouling is diagnosed when one of said evaporator conductance-area product and 
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IS 

3 said evaporator heat exchanger approach exceeds said predetermined threshold value 

4 for said CQ. 

1 32. The fault detection system according to claim 27 wherein condenser 

2 tube fouling is diagnosed when one of said condenser conductance-area product and 

3 said condenser heat exchanger approach exceeds said predetermined threshold value 

4 for said CQ. 

1 33. The fault detection system according to claim 27 wherein an internal 

2 compressor internal fault is diagnosed when one of said compressor coefficient of 

3 performance, said isentropic efficiency, and compressor power draw exceeds said 

4 predetermined threshold value for said CQ. 

1 34. The fault detection system according to claim 27 wherein a motor 

2 transmission fault is diagnosed when one of said overall coefficient of performance, 

3 said motor/transmission efficiency, and motor power draw exceeds said 

4 predetermined threshold value for said CQ. 

1 35. The method according to claim 27 wherein said predetermined 

2 threshold includes upper and lower critical value levels for each said characteristic 

3 quantity, said upper and lower critical value levels being determined in accordance 
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with a rated precision error of a sensor used to acquire a selected said measured input 
used to calculate said characteristic quantities. 

36. The fault detection system according to claim 26 further comprising: 
means for producing interpolated base-case CQ values from base-case chiller 
for a given set of measured inputs; 

wherein said fault detecting means detects a fault when a difference between 
said interpolated base-case CQ values and CQ values for other ones of the chillers 
exceeds a predetermined threshold range for at least one of said CQ values. 
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FIG. 2 




8 Measured 
Inputs 



FIG. 3 

Procedure LMTDe_(T_chws 9 T_chwr 5 T__evap:LMTDe 5 T__we 5 ln_arg_e) 
T_we = (T__chws +T_chwr)/2 

ln_arg_e = max(0,(T_chwr-T__evap)/(T__chws-T_evap)) 
If ((T___chws = T_chwr) OR (ln__arg_e=0)) Then 
LMTDe = T_we - Tevap 
Else 

If (T evap > T_chws) Then Call Error 
If (T_cfaws > T_chwr) Then Call Error 

LMTDereal = ((T_chwr-T _evap)-(T chws-T evap))/lB(ln arg 
LMTDe = max(0,LMTDe_real) 

Endlf 
End 



FIG. 4 

Procedure LMTD_c(T_cws 5 T_cwr,T_cond:LMTDc ? T_wc,ln_arg_c) 
T_wc = (T_cws + T_cwr)/2 

!n_arg_c = max(0 3 (T_cond-T_cws)/(T__cond-T___cwr)) 
If (T_cwr = T__cws) OR (ln_arg_c =0) Then 
LMTDc = T_cond - T_wc 
Else 

If (T__cond < T_cwr) Then Call Error 
If (T__cwr < T_cws) Then Call Error 

LMTDc_real - ((T_cond-T_cws)~(T__cond-T_cwr)) / ln(ln_arg_c) 
LMTDc = max(0,LMTDc_real) 

Endlf 
End 
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time = tableran# 

T__cws = lookup(time, #T_cws) n F" 
Gpmcw = lookup(time, #Gpmcw) ?5 gpm" 
Tcond = lookup(time, #T_cond) f, F" 



FIG. 6 

T_chws = lookup(time, #T_chws) ?, F" 
T_chwr = lookup(time 3 #T_cfawr)"F n 
Gpmchw = Iookup(tirae, #Gpmchw) "gpm" 
T_evap = lookup(time, #T_evap) ,f F" 



FIG. 7 

T_2 = lookup(time 3 #T_2) "F" 
Power = lookup(time 3 #Power) "kW f ' 



FIG. 8 

Q cvap = Gpmchw * ft3hr\gpm * dens chw * Cpw * (T chwr - T chws)/ 12000 "ton" 
Q evap - UAe * LMTDe ,{12000 "ton" 

Qevap = GprnJR. * ft3hr\gpm * dens_R* (h_l - h_4)/12000 "toe" 

h_l = Entha!py(R22, x = 1, P = P_evap) "Btu/lbm" 

v_l = Volume(R22, x = 1, T = T__evap) "ft A 3/lbm" 

T_evap = Temperature(R22,x = 0,P = P_cvap) "F" 

APPRe = T_chws - T evap "F" 

CH WTD = T chwr - T chws "F" 



FIG. 9 

Pcomp - GpmR * ft3hr\gpm * dens R * (h_2 - h i)/ 12000 * kW\ton "kW" 
P_comp = Pideal/eta "kW" 

P ideal = Gpm R * ft3hr\gpm * dens_R * (h_2s - h_l)/ 12000 *kW\ton "kW" 
s_l = Entropy(R22, x = 1, P = P_evap) "Btu/lbm-R" 
s_2s = s_l "Btu/lbm-R" 

h_2s - Enthalpy(R22, s = s_2s, P = P_cond) "Btu/Ibm" 
COP_overall = Q_evap /(Power / kW\ton) 
COP_comp = Q_evap /(Pcomp / kWXton) 
COP_ideal = Q evap /(P_ideal / kWNton) 
EFFMOT = P_comp/Powcr 
P_ratio = P_cond/P_evap 



FIG. 10 

"CONDENSER EQUATIONS:" 

Q_cond = Gprncw * ft3for\gpnn * dens_cw * Cpw * (T_cwr - T ews)/ 12000 "ton" 
Qcond - UAc * LMTDc /1 2000 "ton" 

Q cond = Gpm R * ft3hr\gpm * dens R * (h_2 - h 3)/ 12000 "ton" 

h_3 = Enthalpy(R22 3 x = 0, P = P_cond) "Btu/lbm" 

v_3 = Volume(R22, x = 0, P = P_cond) "ft A 3/lbm" 

T cond = Temperature(R22, x = 0, P = P cond) "F" 

T_2 = Temperature(R22, h = h_2, P = P_cond) "F" 

v_2 = Volume(R22, h = h_2, P = P_cond) "ft A 3/lbm" 

APPRc - T cond - T cwr "F" 

CWTD = Tcwr - T ews "F" 

FIG. 1 I 

"EXPANSION VALVE EQUATION:" 
h 4 = h 3 "Btu/lbm" 



FIG. 12 

"PROPERTY & CONVERSION DATA:" 

dcns_chw = Density( Water, T=(T_chws+T_chwr)/2, P=14.7) "lbm/ft A 3" 

denscw = Density(Water, T=(T_cws+T_cwr)/2, P=14.7) "lbm/ft A 3" 

dens_R = Density(R22, T=(T_cond+T_evap)/2, P=(P_cond +P_evap)/2) "lbm/ft A 3" 

Cpw = 1.0 "Btu/lbm-R" 

ft3hr\gpm = ccnvert(gal,ft3)/convert(min,hr) 

kWYton = convert(ton,kW) 
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FIG. 15 
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FIG. 20 
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FIG. 21 
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FIG. 22 
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FIG. 24 
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FIG. 25 
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FIG. 26 
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